Polarization-selective excitation of N-V centers in diamond 
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The nitrogen-vacancy (N-V) center in diamond is promising as an electron spin qubit due to its 
long-lived coherence and optical addressability. The ground state is a spin triplet with two levels 
{rUs = ±1) degenerate at zero magnetic field. Polarization-selective microwave excitation is an 
attractive method to address the spin transitions independently, since this allows operation down to 
zero magnetic field. Using a resonator designed to produce circularly polarized microwaves, we have 
investigated the polarization selection rules of the N-V center. We first apply this technique to N-V 
ensembles in [100] and [lll]-oriented samples. Next, we demonstrate an imaging technique, based 
on optical polarization dependence, that allows rapid identification of the orientations of many single 
N-V centers. Finally, we test the microwave polarization selection rules of individual N-V centers 
of known orientation. 



I. INTRODUCTION 

Nitrogen-vacancy (N-V) centers in diamond are a 
promising system for implementing and testing quantum 
information processing protocols in solids. Spin states 
of single centers can be initialized and detected opti- 
calljsii^i^i^i^, and this capability has been extended to 
controlled coupling between N-V centers and other spins 
(^^C,7V)Siii^i^ii^. A special feature of this system is the 
spin-triplet structure of the ground states, with individ- 
ual transitions that can in theory be manipulated inde- 
pendently using microwave polarization selection rules. 
The ground state of the negatively charged N-V center 
is known to have a {^A) spin triplet structurei^ which, 
neglecting hyperfine interactions, is described by the fol- 
lowing electron spin Hamiltoniani^,: 
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where B is the magnetic field, g/3 — 27.98 GHz/T and 
D = 2.88 GHz is the zero- field splitting between the 
lowest-energy nis = sublevel, and the nis = ±1 sub- 
levels which are degenerate at zero magnetic field. In 
typical optically-detected magnetic resonance (ODMR) 
experiments, the microwave transitions are driven by lin- 
early polarized fields, and individual transitions are se- 
lected by applying a magnetic field to lift the to^ = ±1 
degeneracy. To avoid spin mixing, the magnetic field is 
usually applied along the quantization axis, which can be 
along any of the four (111) crystal axes. 

An alternative method to drive the spin transitions 
selectively is to take advantage of polarization selection 
rules, since the transitions from = to = ±1 
occur with circularly polarized microwave excitation. 
Using polarization, rather than frequency, selective ex- 
citation at zero magnetic field becomes possible, and 
at weak magnetic fields where transitions still overlap 
{e.g. B < 5G) the selectivity can be improved. At 
small magnetic fields and under weak strain or electric 
field, it is expected that simultaneous cycling and non- 



spin-conserving optical transitions can be realized, useful 
for optical readout and manipulation, respectivel y "'^'^d'^ . 
Therefore, high selectivity and weak magnetic fields are 
highly desired for improving the prospects of N-V-based 
quantum information processing. 

Here we demonstrate selective excitation of ensembles 
and individual N-V centers using circularly polarized mi- 
crowaves for the spin transitions, and the laser polariza- 
tion with respect to the N-V axis for the optical transi- 
tions. 



II. EXPERIMENTAL SETUP 

To test the selection rules experimentally, we per- 
formed optically-detected magnetic resonance (ODMR) 
experiments using a microstrip resonator designed to 
produce circularly polarized microwaves. The cavity is 
formed by two crossed microstrip resonators with or- 
thogonal modes which generate microwaves with any po- 
larization (from linear to circular) when the two inputs 
are driven with a particular phase difference (0° to 90°). 
Our cavity has a full-width half maximum (FWHM) res- 
onance of 30 MHz with a quality factor oi Q — 100. To 
avoid distortions from the resonance shape of the cav- 
ity, the measurements were kept within the FWHM re- 
sponse of the cavity. After positioning the sample above 
the cavity center, the resonance frequency of the cav- 
ity was trimmed by placing a low-loss dielectric above 
it. Three different microwave polarization configurations 
were used: linearly polarized, clockwise (cr~) and coun- 
terclockwise (cr"*") circular polarizations with respect to 
the static magnetic field. The static magnetic field was 
produced by a permanent magnet placed in a linear 
translator close to the sample. The magnetic field di- 
rection was set along the optical detection axis (perpen- 
dicular to sample surface). Apart from the microwave 
polarization control, the techniques used here are simi- 
lar to those used in previous room-temperature ODMR 
experimentsiiii. Fig. [H^a) shows a schematic for our ex- 
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FIG. 1: ODMR experiment on (100) sample a. 
Schematic drawing sliowing tlie experimental positioning of 
the sample relative to the cavity and the magnetic field direc- 
tion, which is also the rotation axis for the ac magnetic field, 
b. ODMR ratio (ct^/ct" amplitude) vs. microwave power in- 
tensity (points) , fitted empirical equation as described in text 
(solid line), and saturation power (dashed line), c. Ensemble 
ODMR spectra for three microwave polarizations: g~ (blue 
circles), cr^ (red diamonds) and linear (green squares). The 
solid lines are fits of gaussian peaks with adjustable ampli- 
tudes and 5 MHz width. The inset is a schematic drawing 
showing the 4 possible N-V orientations. 



perimental setup, showing the sample position above the 
cavity, as well as the cavity design. We used a confocal 
setup in which an excitation laser (532 nm) was focused 
to a sub-micron spot inside the sample for non-resonant 
excitation of the N-V centers. This setup also allowed for 
the control of the laser polarization incident on the sam- 
ple. Laser excitation produces a spin polarization, pref- 
erentially populating the rus = ground state. When the 
applied microwave field with the appropriate polarization 
is resonant with one of the spin transitions, the spin po- 
larization decreases, producing a measurable decrease in 
the photoluminescence intensity. 



III. RESULTS AND DISCUSSION 

A. [100] Oriented Sample: Ensemble Measurements 

The first set of measurements (Fig. [T|) was performed 
on a [100] oriented sample with a relatively high (~ 
10^ /^m~^) concentration of N-V centers near the surface. 
For this crystal orientation, none of the four possible N- 

V orientations (Fig. [U^c), inset) have their quantization 
axes parallel to the surface normal, and hence a circular 
polarization in the laboratory frame of reference has an 
elliptical projection onto a plane perpendicular to the N- 

V axis. However the ODMR spectrum is the simplest for 
this crystal orientation because for each N-V orientation 
the absolute angle between the quantization axis and the 
surface normal is the same. The normalized transition 
probabilities between mg = to uis ± 1 are calculated 
theoretically to be (2 + \/3)/12 and (2 - \/3)/12, respec- 
tively. 

For a quantitative verification of the selection rules we 
would ideally perform measurements in a regime where 
the microwave power dependence of the ODMR spec- 
trum is linear. The ODMR signal dependence on the 
microwave power excitation can be phenomenologically 
described hy A = A^^xli'^ + Psat/P), where Psat is the 
-3dB saturation level. Figure [TJb) shows a measurement 
of the resonance amplitude ratio between the two tran- 
sitions {vris = — > vrLs = -f 1 and rn^ = — > = — 1) 
as a function of applied microwave power. The solid 
line is a fit using the ratio between the two transitions, 
{P + P^J^)/{P + P^^ ) (for the cr± microwave po- 
larization), where in the low power limit the ratio is 

PSl/PS,t = (2-hV3)/(2- V3). With these assumptions 
the saturation power for ct^ is ^ — 49dBm. Our measure- 
ments were done at — 40dBm, since the time required to 
perform the experiment below the saturation limit at the 
same signal to noise ratio would be prohibitively long. 

Figure [U^c) shows ESR spectra measured for three mi- 
crowave polarizations, showing two peaks corresponding 
to the rus = ±1 transitions. The solid lines are the cal- 
culated spectra solving Eq. [T] over the 4 possible N-V 
orientations, considering an inhomogeneous broadening 
of 7.5 MHz (obtained from the measured data) for each 
transition. The model does not include the saturation 
behavior. Circularly polarized microwave excitation pro- 
duced a constrast ratio greater than 75% between the 
two ODMR peaks. In the linearly polarized spectrum 
(shifted for clarity), the slight asymmetry is due to the 
cavity response. 



B. [Ill] Oriented Sample: Ensemble Measurements 

Next, we describe measurements performed on a [111]- 
oriented sample with a much lower N-V concentration, 
below l^m~'^. For this crystal orientation, from four 
possible N-V orientations (see Fig. (Ha)), one of them 
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FIG. 2: Ensemble ODMR spectra for [111] oriented 
sample showing hyperfine interaction, a. Schematic 
drawing showing the relation between the 4 possible N-V ori- 
entations and the magnetic field direction, b. Color plot 
showing ODMR spectra vs. magnetic field (horizontal axis, 
from 1 to 4 gauss) and microwave frequency (vertical axis, 
50MHz range) for linear, a~ and a'^ microwave polarizations, 
comparing data (left) and simulation (right) as described in 
the text. The simulation uses an inhomogeneous broadening 
of 1.2MHz to fit the data. The color scale indicates the peak 
depth, in percent, relative to the off-resonant case. The aver- 
age number of counts, for the three measurements, is 4 x 10^. 



([Ill], labeled "NVl") is normal to the sample surface 
and is expected theoretically to exhibit perfect selection 
rules. For the other orientations the projected microwave 
field is elliptical in the N-V reference frame, and the cir- 
cular selection rules are not fully observed. Thus, this 
measurement can reveal the orientation of a particular 
center. For this experimental configuration, the Zeeman 
splitting depends on the N-V orientation, decreasing by 
a factor of approximately 3 for N-V orientations other 
than [111]. 

Due to the high purity of this sample, hyperfine struc- 
ture can be resolved in the ODMR spectrum. For a sub- 
stitutional ^'^N with nuclear spin 7=1, the Hamilto- 
nian terms describing the nuclear Zeeman, hyperfine and 
quadrupole interactions arei^: 
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Ht = P{ll - - (.g„/?„)B • I + ... 

+A\\Sj,, + \a^{s-i+ + s+r).. 



(2) 



where /3„ is the nuclear magneton, A\\ and A±^ are the 
hyperfine constants and P is the quadrupole splitting. 
Previous worki^ii^iil has established these constants to 
be: A\\ = 2.3 MHz, Ax_ = 2.1MHz, 5„/3„ = 3.07MHz/T 
and P = —5.04 MHz. For small static magnetic fields, 
the hyperfine interaction splits each of the two ODMR 
peaks into a triplet with splittings of approximately A|| = 



2.3 MHz. In this limit, each component of a triplet follows 
approximately the original circular polarization selection 
rules. 

The ensemble ODMR spectra in Fig. ^h) , which show 
the hyperfine structure, were obtained by removing the 
pinhole from the collection optics to reduce spatial se- 
lectivity, and exciting with small optical and microwave 
powers. Measured and calculated spectra using Eq. [2] 
are shown. The calculation includes all possible tran- 
sitions, weighted by their predicted strengths, using a 
lorcntz function (1.2 MHz linewidth) to account for in- 
homogeneous broadening. For each panel in [D^b) , the 
color scale represents the ESR intensity as a function 
of the static magnetic field (x-axis) and the microwave 
frequency (y-axis) for the three different microwave po- 
larizations. Two different slopes are observed for the 
Zeeman splitting, corresponding to N-V centers oriented 
along [111] (steeper slope) and the other three orienta- 
tions. These results exhibit good agreement between the- 
ory and experiment. For the N-V centers oriented along 
[111], the measured contrast ratio for vs. a~ excita- 
tion is approximately 9. 



C. Optical determination of defect orientation 

For work on single N-V centers, it is helpful to have a 
simple method to determine the orientation of the N-V 
axis. Here, we demonstrate a method based on opti- 
cal polarization dependence. Since the dipole moments 
of the optical transitions are perpendicular to the N-V 
axis, the optical polarization anisotropy depends on ori- 
entationi^. For a [111] sample, all four orientations can 
be distinguished optically, allowing rapid identification 
of many N-V centers through polarized scanning confo- 
cal imaging. For effecient representation of the data, we 
combine many images obtained with different excitation 
polarizations into a single image by encoding polarization 
into color according to the following filter function: 



/„(0) = a + sin 2(0 4- (/lo + m:/3) , 



(3) 



where /„ is the multiplication factor for color channel 
n = 0,1,2, (j) is the polarization angle, a is an offset 
chosen as small as possible to maximize color saturation 
without producing negative pixel values, and 0o is an 
angle offset which controls hue. This procedure is il- 
lustrated in Fig. EI a) , where the first three panels show 
images obtained at the particular polarizations for which 
N-V orientations NV2, NV3, and NV4 have photolumi- 
nescence minima. The colors of these individual images 
correspond to the encoding in Eq. [3l The insets show 
the relative orientation of the laser polarization (i) with 
respect to these three N-V axes. The minima occur when 
the projection of the N-V axis onto the sample surface 
is parallel to the electric field of the optical excitation. 
The last panel shows the combined image in which the 
four orientations can be clearly distinguished by color. 
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FIG. 3: Determination of N-V orientation through optical polarization dependence, a. Scanning confocal photo- 
luminescence images showing single N-V centers, with laser polarization encoded into the three color channels MYC (magenta, 
yellow, and cyan). The first three images were measured for particular polarizations in which orientations designated as NV2, 
NV3, and NV4 have photoluminescence minima, with the electric field nearly parallel to the N-V axis. A representative N-V 
center of each orientation is also labeled. The fourth image is a sum over many measurement polarization angles, in which orien- 
tations NVl, NV2, NV3 and NV4 appear as white, green, blue and red, respectively, b. Luminescence intensity vs. excitation 
polarization for selected single N-V centers with orientations NVl (black triangles, oriented normal to sample surface), NV2 
(green squares), NV3 (blue diamonds), and NV4 (red circles). These measurements were performed at weak excitation power 
to avoid saturation of the optical transition, which distorts the polarization dependence, c. Excitation power dependence for 
a single N-V center with NVl orientation, fitted to a simple model (solid line) . 



Based on photon correlation measurements, we can un- 
ambiguously associate the individual spots with single 
N-V centers. 

The optical polarization dependence of four represen- 
tative N-V centers is shown in Fig[3]Jb), fitted to a simple 
model. This model includes the effect of N-V orientation 
on both the excitation rate and the photoluminescence 
collection efficiency. For an N-V center oriented normal 
to the sample surface (NVl), the two dipole moments 
(perpendicular to the N-V axis) can be excited equally 
by the laser. Experimentally, for this orientation we ob- 
serve that the photoluminescence is approximately unpo- 
larized, independent of the excitation polarization (not 
shown). Therefore in the model we assume that the two 
dipole moments of the N-V center radiate with equal in- 
tensity, independent of which dipole moment is excited 
by the laser, a situation that would occur if relaxation 
between the excited states is fast at room temperature. 



When both polarizations are detected as in Fig[3fb), the 
other orientations (NV2, NV3, NV4) produce a weaker 
photoluminescence signal even at maximum excitation ef- 
ficiency since the collection efficiency is reduced for this 
geometry. The model also includes saturation of the op- 
tical transitions (measured in Fig[3fc)) and asymmetry 
in the collection efficiency in our setup for light polar- 
ized parallel or perpendicular to the laser polarization. 
Taking all of these factors into account, we obtain good 
agreement between theory and experiment. 



D. [Ill] Oriented Sample: Single NV 
Measurements 

Finally, using the method above to locate single N-V 
centers of known orientation, we proceed to test their mi- 
crowave polarization selection rules. The ODMR spectra 
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FIG. 4: Single N-V ODMR measurements, a. Signal for 
a single N-V center oriented along [111] (NVl) for (blue 
circles), ct"*" (red diamonds), and linear (green squares) mi- 
crowave polarizations. The solid lines are fits using a sum of 
two gaussians with adjustable amplitude, b. Similar mea- 
surements performed on a single center with a different ori- 
entation (NV2), fitted using a single gaussian. The Zeeman 
splitting is reduced in this case because of the smaller projec- 
tion (cos 6 — 1/3) of the magnetic field onto the N-V axis. 



for single N-V centers with NVl and NV2 orientations 
are shown in Fig. [H with data (filled symbols) and fits 
(solid lines). In Fig. [4l^a) the fit used a superposition of 
two gaussian curves corresponding to the two microwave 
transitions, while in Fig. [DJb) the fit used a single gaus- 
sian since, for the smaller Zeeman splitting, the transi- 
tions are unresolved. For NVl we expect perfect selectiv- 
ity, characterized by a forbidden transition from m = 
to m = —1 under cr"*" excitation. Experimentally, we 
observe a residual peak in this case with approximately 
15% relative intensity which is most likely due to imper- 
fect generation of circularly polarized microwave fields 



above the cavity, or even the parasitic effect of the trim- 
ming dielectric, sample shape and microscope objective. 
Simulations predict that if the measurement location is 
just a few millimeters away from the cavity center, some 
ellipticity can occur. For the NV2 orientation, although 
the Zeeman splitting is not resolved, exciting with cir- 
cular microwave polarization causes a narrowing of the 
transition linewidth compared with linearly polarized ex- 
citation, and a shift of the peak center is also observed, 
showing that the circular polarization favors one transi- 
tion. Theoretically, we expect the Zeeman splitting for 
the NV2 orientation to be smaller by a factor of 3. We 
also expect reduced polarization selectivity due to the 
large angle between the N-V axis and the microwave po- 
larization axis, with an additional contribution from the 
cross-polarized transition of 1/4 amplitude. 



IV. CONCLUSIONS 

In summary, we have experimentally verified spin po- 
larization selection rules for N-V centers in diamond with 
various crystal and N-V orientations using circularly po- 
larized microwave excitation. The results, including hy- 
perfine structure, are consistent with theory. We have 
also demonstrated an efficient optical method for deter- 
mining orientations of single N-V centers, and demon- 
strated microwave polarization selection rules in single 
N-V centers of known orientation. By optimizing the rel- 
ative amplitudes and phases of the signals used to drive 
the two inputs of the microwave cavity, a much higher 
selectivity is possible. Even for the case of the [100] ori- 
ented sample, for a single N-V center, elliptically polar- 
ized excitation could in principle give perfect selectively. 
Therefore we expect that microwave polarization selec- 
tion rules will prove quite useful for spin manipulation, 
initialization and tomography of single N-V centers when 
improved selectively or operation at zero static magnetic 
field are needed. 
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